Mon. Not. R. Astron. Soc. 0Q0,rrifl3l(2Q08) Printed 9 May 2008 (MN LSTbX style file v2.2) 



The ionizing radiation from massive stars and its impact on HII 
regions: results from modern model atmospheres 



S. Simon-Diaz^'^ * and G. Stasinska^ 

' ^LUTH, Observatoire de Paris, CNRS, Universite Paris Diderot; 5 Place Jules Janssen, 92190 Meudon, France 
^Geneva Observatory, University of Geneva, 51 chemin des Maillettes, 1290 Sauverny, Switzerland 



Accepted//Received 



ABSTRACT 

We present a detailed comparison of the ionizing spectral energy distributions (SEDs) pre- 
dicted by four modern stellar atmosphere codes, tlusty, cmfgen, WM-basic, and fastwind. We 
consider three sets of stellar parameters representing a late O-type dwarf (09.5 V), a mid O- 
type (07 V) dwarf, and an early O-type dwarf (05.5 V). We explore two different possibilities 
for such a comparison, following what we called evolutionary and observational approaches: 
in the evolutionary approach one compares the SEDs of stars defined by the same values of 
Tgff and \ogg; in the observational approach the models to be compared do not necessarily 
have the same T^ff and log g, but produce similar H and He i-n optical lines. We find that there 
is a better agreement, in terms of Q(li the ratio Q(lie "VQCH and the shape of the SEDs 
predicted by the four codes in the spectral range between 13 and 30 eV, when models are 
compared following the observational approach. However, even in this case, large differences 
are found at higher energies. 

We then discuss how the differences in the SEDs may affect the overall properties of sur- 
rounding nebulae in terms of temperature and ionization structure. We find that the effect over 
the nebular temperature is not larger than 300-350 K. Contrarily, the different SEDs produce 
significantly different nebular ionization structures. This will lead to important consequences 
on the establishment of the ionization correction factors that are used in the abundance deter- 
mination of H II regions, as well as in the characterization of the ionizing stellar population 
from nebular line ratios. 

Key words: Stars: early-type - Stars: fundamental parameters - Stars: atmospheres - ISM: 
H II regions 



1 INTRODUCTION 

The intense far ultraviolet (FUV) radiation emitted by early OB- 
type stars ionizes the interstellar medium, generating the so-called 
Hii regions. These regions can be used to derive properties of the 
associated stellar population (e.g. initial mass function, star form- 
ing rate, age), and other properties of the galactic region where 
these are located, such as chemical composition. However, since 
the properties of H ii regions crucially depend on the ionizing spec- 
tral energy distribution (SED) of the massive star population, and 
this part of the stellar flux is generally unaccesible to direct ob- 
servations, the predictions resulting from massive star atmosphere 
codes are a crucial ingredient. 

The outer layers of blue luminous stars are characterized by 
strong non-LTE conditions, spherically extended geometries, and 
the effect of hundreds of thousands of fiux absorbing metal lines 
present in the FUV and UV spectral ranges (producing the so- 
called line-blanketing and line blocking effects, as well as the de- 
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velopment of radiatively driven stellar winds). All the above ef- 
fects must be taken into account when modeling the atmospheres 
of these stars. An ideal model atmosphere code should consider all 
of these problems in a detailed manner; however, this would require 
an enormous computational effort. To avoid this, depending on the 
specific study one wants to treat, some of the physical processes 
can be relaxed without affecting the reliability of the model output. 

In the last decades, a great effort has been devoted to the devel- 
opment of stellar atmosphere codes for hot massive stars. The ad- 
vent of the new generation of non-LTE, line blanketed mode l atmo - 
sphere codes, either plane-parallel (tlusty, Hubeny & L anz. | l995h . 
or spherically expanded (fastwind, Santolaya-Rey et al. 119971, Puis 
et al. 200 5; cmfgen, Hillier & Miller 1998; WM-basic, Pauldrach et 
al. l200ll) is already a fact. While cmfgen aims to be the most ex- 
act one (and hence, the most time consuming), the other three use 
specific approximations in the calculation of the stellar atmosphere 
structure and of the SED, reducing then the computational time. 

New generation stellar atmosphere models, including a more 
realistic description of the physical processes characterizing the 
stellar atmosphere, produce quite different ionizing SEDs than the 
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Table 1. Summary of stellar atmosphere models used in this work (C. W, F, 
and T labeling cmfgen, wu-basic, fastwind, and tlusty, respectively). 



Models 




logg 


R 


logM 






(kK) 


(dex) 


{Re) 


(Moyr-') 


(kms-') 


Cl.Wl.Fl.Tlt 


30.0 


4.0 


7.0 


-7.28 


2000. 


C2, W2, F2, 72+ 


35.0 


4.0 


9.0 


-7.12 


2000. 


C3, W3, F3, T3+ 


40.0 


4.0 


10.0 


-6.93 


2400. 


F4 


31.0 


4.1 


7.0 


-7.28 


2000. 


W4 


31.0 


4.0 


7.0 


-7.28 


2000. 


T4 


30.5 


4.1 








F5 


36.0 


3.9 


9.0 


-7.12 


2000. 


W5 


36.0 


4.0 


9.0 


-7.12 


2000. 


T5 


35.5 


4.1 








F6 


41.0 


4.0 


10.0 


-6.93 


2400. 



' R. M, and i>„ only for cmfgen, WM-basic, and fastwind models. 

previ ous plane-parallel, hydrostatic mode ls (eith er LTE by Kurucz 
I l99ll or non -LTE by e.g. Mihalas & Auer ll97(]|. Kunze et al. ll992l. 
Kunze Il994h . Some notes on this, and on the consequences on 
the ioniza tion structure o f H ii regions, can be fo und in Gabler et 
al. ( (l989j) , Naja rro et al. ( Il996ll . Sellmai er et a l. ( Il996h . R ubin e t 
al. 



1991 l2007t) , Stasinsk a & S chaerer jl997h . Schaerer ( I2OO0I) , 
Martin-Hernandez et al. Although the new p redictions 

seem s to go in the right direction (viz. Giveon et al. 2002, Morisset 
et al. l2004h non-negligible differences can still be fou nd between 
th e vario us stellar codes (see e.g. Mokiem et al. |2004| . Martins et 
al. bOOl Puis etal. l2005h . 

In this paper, we explore in detail the differences between the 
predictions in the FUV from four modem stellar atmosphere codes 
(tlusty, cmfgen, WM-basic and fastwind), and study their effects 
on H n regions spectra. In Sect. |2l we present a set of stellar at- 
mospheres constructed specifically for this purpose, and show that 
two different approaches can be followed to compare the SEDs: the 
evolutionary approach, in which the same values of T^f{ and log g 
are used for models computed with the various stellar atmosphere 
codes; and the observational approach, in which the models to be 
compared do not necessarily have the same T^fj and log g, but pro- 
duce similar H and He i-ii optical lines. In Sect. [3] by using sim- 
ple, ab initio, H n region models, we study how the differences in 
the SEDs may affect the overall properties of surrounding nebulae 
in terms of temperature, ionization structure and other astrophys- 
ical diagnostics. Finally, Sect. |4] summarizes the main results and 
presents some prospects. 



2 COMPARISON OF IONIZING SEDS FROM VARIOUS 
STELLAR ATMOSPHERE CODES 

For the reader unfamiliar with stellar atmosphere modelling, the 
main characteristics of the codes tlusty, cmfgen, WM-basic and 
FASTWIND are presented in the appendix (accessible on-line only). 



dex and effective temperatures T^ff = 30000, 35000, and 40000 kQ, 
and computed cmfgen, fastwind and WM-basic models for each pair 
Teif-logg. In order to exclude from our study the possibility that 
discrepancies found between the stellar atmosphere codes calcula- 
tions could be due to differences in the considered stellar and wind 
parameters and/or stellar abundances, those models were calculated 
taking care of using exactly the same set of parameters and abun- 
dances. 

As part of the initial set of models, three tlusty models, taken 
from the ostar2002 gric[|by Lanz & Hubeny l l2003l) . were also con- 
sidered (those labeled as G30000g400vl0, G35000g400vI0, and 
G40000g400vI0, respectively). 

Finally, another set of fastwind, wu-basic and tlusty mod- 
els, with slightly modified values for the effective temperatures and 
gravities, was also computed. This was necessary for the approach 
presented in Sect. l2.4l 

In all the cases, the metallicity was considered to be sol ar (fol - 
lowing the set of abundances derived by Grevesse & Sauval. ll998h . 
and a microturbulence =10kms"' was assumed in the model cal- 
culations. 



2.2 Methodology 

Two different approaches were followed to compare the SEDs from 
the various stellar atmosphere codes. Basically, these are charac- 
terized by the methodology used to establish the stellar parameters 
(Teff , log g, and log L). We call them evolutionary and observational 
approaches. 

In the evolutionary approach the stellar parameters character- 
izing a star are fixed by the output of stellar evolutionary codes. 
This approach is mainly used b y stell ar population synthe sis cod es, 
such as STARBURST99 (Leitherer l 19991 see also Smith et al. I2OO2I) or 
PEGASE (Fioc & Rocca-Volmerange ll999l) . In this case, for a given 
initial stellar mass and evolutionary time, a star is defined by the 
Teif, \ogg, and logL predicted by the stellar evolu tion codes (e.g. 
Schaller et al. 1992 ; Meynet & Maeder 2000 ll2005l) . 

On the other hand, when spectroscopic and photometric 
observations of the star are available, the stellar parameters can be 
derived by means of what we call an observational approach. In 
this case, the stellar and wind parameters are determined by fitting 
appropriate stellar line features using stellar atmosphere models. 
This is a long established way to analyze O-type star spectra. The 
optical H and He i-ii lines have been mainly used to determine the 
stellar T^ff and logg (viz. Herrero et al. 1992, 2002; Repolust et al. 
2004), though the UV spectra of these stars also provide valuable 
information about the physical properties of their winds, such 
as terminal wind velocities and mass loss rates (see Kudritzki & 
Puis 200?). In addition, in the observational approach, the stellar 
luminosity is obtained from the stellar absolute visual magnitude 
(M,). 

For a quantitative comparison of ionizing SEDs, we used the 
quantity Qe, indicating the number of photons with energy greater 
than E=hy. If the emergent stellar flux is given by Hy(A), these 
quantity is defined as follows: 



2.1 Choice of stellar atmosphere models 

Table [T] summarizes the stellar atmosphere models used for our 
study. Initially, we considered three stars with a gravity logg = 4.0 



' Following the spectral type - T^tf calibration by Martins et al. j2005l) . 
these models would approximately represent a late O-type (09.5 V), a mid 
O-type (07 V), and an early O-type (05.5 V) dwarf, respectively. 
^ See also |http://nova.astro.umd.edu/| 
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Figure 1. (Upper panels) Spectral energy distiibution in the range 10-60 eV as predicted by the various stellar atmosphere codes, cmfgen, tlusty, v/M-hasic, 
and FASTWiND predictions are plotted as black, red, blue, and green curves, respectively. The fluxes from WM-hask, cmfgen and tlusty have been re-mapped 
to the FASTWE^D frequency grid for a better comparison. (Lower panels) Comparison of the number Qe of photons with energy greater than E relative to the 
con'esponding number for a blackbody with T=T^fi. The location of the ionization edges of various ionic species of H, He, N, O, Ne, S, CI, and Ar with 
ionization energies below 54.4 eV (He + ionization energy) are also shown. 



Qe-1 ^dy' = l6n^R-l ^ di' (1) 

where Ly the emergent stellar luminosity distribution, and R the 
stellar radius. 

2.3 Comparison of SEDs in the evolutionary approach 

In this section, we compare the SEDs resulting from cmfgen, 
TLUSTY, wu-basic, and fastwind when models with exactly the 
same stellar parameters are considered (Tcs, log g and L). Through 
this comparison we can test the implications of using SEDs from 
the various stellar atmosphere codes as input for stellar population 
synthesis codes in order to compute photoionization models of 
giant Hii regions. 

Upper panels in Figure[T]show the corresponding SEDs in the 
range 10-60eV. This qualitative comparison already shows that 
the stellar atmosphere codes produce somewhat different ionizing 
SEDs. To quantify these differences, we first consider the total 
number of H " ionizing photons, and then compare the shape of 
the SEDs for energies E > 13.6 eV. 

2.3.1 Ionizing luminosities 

Table |2] summarizes the number of H" ionizing photons, 
Q(ll°) = Qi3.6, resulting from the various stellar atmosphere mod- 
els. Note that in this case, the stellar radii (used in the calculation 
of the ionizing luminosities by means of Eq. [T]! were fixed once 



Table 2. Resulting number of H ionizing photons (in logarithm) for the 
various stellar atmosphere codes when the same T^g and log g are consid- 
ered (evolutionary approach). 





30000 


35000 


40000 




4.0 


4.0 


4.0 


logL 


38.10 


38.55 


39.80 


C1,C2, C3 


47.33 


48.39 


48.99 


T1,T2, T3 


47.43 


48.44 


49.00 


Wl, W2, W3 


46.97 


48.30 


48.93 


F1,F2, F3 


46.98 


48.26 


48.92 



a characteristic stellar luminosity was assumed for each effective 
temperature. 

The values of log QQi") computed with the various codes for 
the same Teff, logg, logL are found to differ significantly. In the 
worst case, i.e. between WM-basic and tlusty for refr=30000, they 
differ by ~ 0.5 dex. The discrepancy tends to be lower for larger 
Teff , but in the hottest star a difference of 0.08 dex (i.e. a factor of 
1.2) can still be found between fastwind and cmfgen models. Note, 
however, that there is always a good agreement (below 0.04 dex) 
between WM-basic and fastwind models on one hand, and cmfgen 
and TLUSTY models on the other. 
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2.3.2 Shape of the ionizing SEDs 

To study how the shape of the predicted SEDs compares in the 
H Lyman continuurr0 we produced the diagrams shown in lower 
panels of Figure[T] There we compare, for each energy E (eV), the 
number Qe of photons with energy greater than E relative to the 
corresponding number for a blackbody with T=T^ff. Note that all 
spectra were scaled to have the same value for 2i3.6- 

These diagrams can also be used to compare the number of 
ionizing photons for a certain ionic species, X ', relative to the total 
number of hydrogen ionizing photons - Q(K ')/2(H ") This is the 
major factor determining the relative populations of different ions 
of a given element in an H ii region, although its effect is modulated 
by the recombination and charge exchange coefficients. 

Since we are interested on the impact that the ionizing radi- 
ation from massive stars has on Hii regions, we only considered 
ions of those elements which lines are commonly observed in 
Hii regions (i.e. H, He, N, O, Ne, S, CI, Ar). The corresponding 
ionization potentials are indicated as vertical lines in Figure[T] Note 
that, for the range of effective temperatures we are considering, all 
the SEDs show an abrupt decrease at energies >54.4 eV. The effect 
of this part of the stellar flux on the ionization of the surrounding 
material in the case of H ii regions is practically negligible, so we 
restrict our study to ions with ionization potential <54.4 eV. 

Various conclusions concerning the shape of the SEDs can be 
emphasized from the inspection of Figure[T] 

30000 K star: In this case fastwind and WM-basic models predict 
a larger H Lyman jump than the other two codes. This is consis- 
tent with the fact that a lower number of H " ionizing photons is 
resulting from fastwind and wu-basic models (see above). 
Good agreement is found for all the codes except tlusty in the 
region between the H" and He" edges. This code is producing a 
somewhat harder flux, and hence larger Q(S and Q(C\ ^) values. 
While CMFGEN and wu-basic predict a similar He " jump, this break 
is larger in tlusty and fastwind SEDs. As a consequence, tlusty 
results in the same Q(lie ") than fastwind and WM-basic, and fast- 
wind in a smaller number of He " ionizing photons (~ 0.2 dex). 
The slope of the SED above the He " edge is larger in tlusty and 
CMFGEN than in wu-basic and fastwind, which results in harder 
fluxes for these later codes for energies E >40 eV. 
It is also interesting to note how the presence of a strong emission 
line in ~ 38 eV (present in WM-basic and cmfgen spectra, but not in 
FASTWIND and tlusty) affects the resulting Q{S and 2(0 ^). 

35000 K star: The agreement between tlusty and cmfgen is quite 
good, with discrepancies not larger than 0. 1 dex in Qe, except in the 
range 30 - 40 eV. This is a region of the stellar flux where metal line 
blocking is important, and differences in the way this is treated (e.g. 
the amount of metal lines included in the blocking calculation), or 
differences in the ionization degree of the elements producing this 
forest of lines can affect the amount of flux emerging from the star 
in this spectral range. The discrepancy found between cmfgen and 
TLUSTY in this spectral range translates in Q(S and 2(0 ^) values 
differing by ~ 0.4 dex. 

Important discrepancies are found between these two codes and 
wu-basic and fastwind in the spectral range where the S^, Cl^, 
He Ar N + ionizing edges are located (~ 23 - 30 eV). Both wm- 



Table 3. Stellar parameters {Ti-g and logg) required to fit fastwind and 
TLUSTY optical H and He lines to those resulting from cmfgen models with 
log g =4.0 and T^^ = 30000, 35000, and 40000 K, respectively. The labels 
indicated in brackets are the same as those used in Table[T] 



CMFGEN 


fastwind 


TLUSTY 


30000,4.0 (CI) 
35000, 4.0 (C2) 
40000, 4.0 (C3) 


31000,4.1 (F4) 
36000, 3.9 (F5) 
41000, 4.0 (F6) 


30500, 4.1 (T4) 
35500, 4.1 (T5) 
40000, 4.0 (T3) 



basic and fastwind result in lower fluxes in this spectral range. For 
example, (2(He '')/2(H ") is ~ 0.3 - 0.4 dex smaller in fastwind and 
wu-basic than in the other codes. In addition, while wu-basic is 
closer to cmfgen and tlusty for larger energies, fastwind is pro- 
ducing a lower flux between 30 and 45 eV. 

In this case WM-basic and fastwind are also resulting is a somewhat 
larger H Lyman edge (again consistent with results in Table|2](. 

40000 K star: Good agreement is found between cmfgen and 
TLUSTY (again except in the range 30-40eV). fastwind predicts a 
somewhat lower flux than these two codes for E>25 eV. On the 
other hand, while wu-basic results in a similar SED as cmfgen for 
E < 35 eV, differences of up to 0.8 dex in Qe are found for larger 
energies. 

2.4 Comparison of SEDs in the observational approacli 

In this section we check if the various stellar atmosphere codes are 
producing similar results for the lines normally used to estimate the 
stellar parameters. Then we compare the SEDs given by models 
which do not have necessarily the same T^s and log g but produce 
similar results for those lines. 

2.4. 1 Predicted stellar HHe optical lines 

Figure |2] shows the comparison of a representative set of H B aimer 
and Hei-n synthetic lines for the various values of (reff, logg), as 
predicted by fastwind, cmfgen, and tlusty. We decided to present 
results in two He i lines to account for the singlet-triplet problem 
(see Najarro et al. 2006, and references therein). WM-basic is not 
included in this comparison since the optical H and He i-ii lines are 
not synthesized by this code. 

While fastwind and cmfgen produce quite similar Hy wings 
(which are the main diagnostic to determine logg), tlusty mod- 
els would need a larger logg to produce the same result. On the 
other hand, while the three codes predict coherent results for the 
He I A 447 1 (triplet) line, this is not the case for the other He i A 4922 
(singlet) line. The caus e of t his discrepancy was already pointed 
out by Najarro et al. who also suggested that the Hei 

triplet lines should be preferred in the spectral analysis due to the 
sensitivity of the Hei singlet transitions to model assumptions. 
Finally, cmfgen and tlusty agree almost perfectly regarding the 
He 11/14541 line, while fastwind models result in a shallower line 
in the T^ff = 30000 and 35000 K cases. Although it is not shown 
here, a similar behavior is also found for the other He ii lines. 

The above-mentioned result implie that the codes are result- 



' Note that shocks were not considered in any of the models and conse- 
quently the emergent stellar flux in the X-ray domain is virtually zero. 



Villamariz & Herrero ilOOd) showed that, except for He ii A 4686, the 
Hen lines are practically insensitive to the adopted value of microturbu- 
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Figure 2. Comparison of some of the optical H, He i, and He u lines commonly used to derive stellar parameters in O-type stars resulting from cmfgen, tlusty, 
and FASTWiND models. Note that only H and He hnes are present in the fastwind synthetic spectra (see appendix). 



ing in a somewhat different temperature and density structure in 
the stellar atmosphere, which in turn affect the strength of the He n 
lines (via the He ionization degree in the line formation region) 
and the wings of the H B aimer lines. On the other hand, the stellar 
parameters that would be derived from the spectroscopic analysis 
of a given observed spectrum are model dependent. To explore the 
differences in the stellar parameters that would be obtained when 
using the various codes we obtained the T^g and logg values re- 
quired by FASTWIND and tlusty to produce similar optical H and 
He i-ii than cmfgen models, taken as reference. These are indicated 
in Tabled 

The discrepancies between the three codes for the 30000 K 
and 35000 K stars are -500- 1000 K in Jeff and 0.1 dex in logg. 
Note that these discrepancies are of the same order as the accuracy 
reached when analyzing observed spectra of OB-type stars with a 
given cod e (~ 1 000 K in T^-g a nd 0. 1 dex in logg — se e e.g. Her- 
rero et al. llOOl Repolust et al. |2004 Simon-Dfaz et al. l2006h . On 
the other hand, the agreement for the 40000 K star case is almost 
perfect. 

In next sections, we take into account these differences for the 
comparison of ionizing SEDs. 

2.4.2 Ionizing luminosities 

Table|4]indicates the number of H " ionizing photons for the models 
summarized in Table[3] The synthetic V magnitudes resulting from 
the models, along with the derived stellar radii, luminosities and 
masses are also indicated. Note that in th is cas e, the stellar radii 
were determined following Kudritzki et al. l ll980h . once an absolute 
visual magnitude (My) was assumed for each set of models. 

In all cases, the agreement between models which results in 
similar H and He i-n lines is better than when the same stellar pa- 
rameters are used to construct the models (see Table[2]l. Note, how- 

lence. Note also that for effective temperatures in the range between 30000 
and 35000 K. the sensitivity of the Hei lines on T^g is quite small. 



Table 4. Number of H ionizing photons for cmfgen, tlusty, fastwind, and 
WM-hasic models with different stellar parameters (T^g and log g), but simi- 
lar H and He optica l lines. The stellar radius was determined following Ku- 
dritzki et al. 1 1980), once an absolute visual magnitude (My) was assumed 
for each set of models. 









Mv = 


-3.5 






Teg, logg 


V 


R/Re 


logL 


M/Mo 


log e(H'') 


CI 


30000, 4.0 


-29.09 


6.3 


38.05 


14.3 


47.28 


T4 


30500, 4.1 


-29.09 


6.2 


38.07 


17.8 


47.45 


F4 


31000, 4.1 


-29.05 


6.3 


38.11 


18.5 


47.40 



My = -4.3 





7'eff, logg 


V 


R/Re 


logL 


M/Ma 


log e(H 


C2 


35000, 4.0 


-29.37 


7.9 


38.52 


22.8 


48.36 


T5 


35500, 4.1 


-29.37 


7.9 


38.54 


28.8 


48.45 


F5 


36000, 3.9 


-29.37 


7.9 


38.57 


18.2 


48.40 



My = -4.8 





Teg, logg 


V 


R/Re 


logL 


M/Ma 


log e(H 


C3 


40000, 4.0 


-29.55 


9.2 


38.88 


31.0 


48.97 


T3 


40000, 4.0 


-29.52 


9.3 


38.89 


31.8 


49.00 


F6 


41000, 4.0 


-29.56 


9.2 


38.92 


30.6 


48.99 



ever, that differences up to 0. 1 dex still remain. These differences 
are partially due to the somewhat different stellar radii (which are a 
consequence of the synthetic V magnitude resulting from the stellar 
models). 

2.4.3 Shape of the ionizing SEDs 

Figure[3]compares the shape of the SEDs for those cmfgen, tlusty 
and FASTWIND models resulting in similar H and He i-n lines. 
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In all the cases, the discrepancies found in the spectral range 
13.6-30 eV when models with the same stellar parameters were 
compared (see Sect. 12. 3| and Fig.[T]( have practically disappeared. 
This is especially noticeable in the 35 and 40 kK star cases. In ad- 
dition, the agreement between the three codes for the H Lyman 
break is quite good. 

Although WM-basic could not be included in the comparison 
presented in Section 12.4.11 (since the optical H and He i-ii lines 
are not sythesized by this code), in view of the abovementioned 
results, we decided to also include the SED of those wu-basic 
models which, having an effective temperature and gravity close to 
those of CMFGEN models, produce similar SEDs as the other codes 
in the range 13.6 - 30 eV. These models correspond to those labeled 
W4, W5, and W3 in Table[Tl 

In general, the situation for higher energies has also im- 
proved; however some important discrepancies still remain. From 
inspection of Fig.[3]it can be remarked that: 

• FASTWiND gives much harder SEDS than the other codes in the 
30000 K case 

• WM-basic results in harder fluxes than the other codes in the 
35000 and 40000 K cases (this discrepancy is especially noticeable 
in the hottest star). 

• TLUSTY and cmfgen SEDs agree quite well, except in the range 
30 - 40 eV. In this energy range, the fastwind SED is closer to that 

of TLUSTY. 

As commented by Puis et al. ( l2005h . the number of metal lines 
between 35 and 50 eV is quite large, and the treatment of the opac- 
ities of the weakest lines affects the SED above 24.6 eV. The differ- 
ent approaches of the line blanketing and wind effects might also 
contribute to the observed discrepancies. In addition, the fact that 
different effective temperatures are being considered in the various 
codes may affect the slope of the SED in the high energy range. 

2.5 The importance of the HHe Une analysis in the 
comparison of SEDs 

In the previous section we have shown that the number of H " ion- 
izing photons, the ratio 2(He'')/2(H''), and the shape of the SEDs 
below 30 eV agree better in the observational approach than in the 
evolutionary approach. This results from the fact that the shape of 
the emergent flux is a consequence of the stellar atmospheric struc- 



ture (which, for a given set of stellar parameters, is influenced by 
the line blanketing and wind-blanketing effects) and the blocking 
by numerous metal lines present in the FUV spectral range. Specifi- 
cally, the size of H " and He " jumps are a function of the population 
of the ground levels of the corresponding ions in the respective con- 
tinuum forming regions. Temperature and electron density mainly 
control the level population of those gro und lev els (though wind 
efl^ects can also contribute, see Gabler et al ll989l , and Najarro et al. 
[199^. For example, in a star with larger eff'ective temperature or 
lower gravity, H and He are expected to be more ionized in the cor- 
responding continuum forming regions, and hence the ground level 
population of H " and He " will be smaller, resulting in smaller H 
and He: jumps. 

On the other hand, the optical H and Hei-ii lines depend 
mainly on the recombination from the corresponding upper ions 
(i.e. H^, He^ and He'^, respectively) which, in turn, is controlled 
by temperature and electron density, again in the region of the stel- 
lar atmosphere where these transitions become optically thin. These 
line formation regions are normally located close to the regions 
where the H" and He" continua are formed, except maybe for H„, 
Hen ^4686, He 1^5875, the strongest lines. 

By fitting the optical H and Hei-n lines, one is fixing the phys- 
ical properties (density and temperature) in the corresponding line 
forming region, and hence in the H " and He " continuum forming 
regions. Consequently, difl"erences in the predicted SEDs, in terms 
of H " Lyman jump and 2(He '')/l2(H ") ratio, and differences in the 
optical H and He ii lines between any couple of models are strongly 
related. 

Therefore, the study of the optical H and Hei-ii stellar lines 
provide very useful additional information to understand possible 
discrepancies in the ionizing stellar SEDs predicted by the various 
stellar atmosphere codes. 

2.6 Use of SpT-reff calibrations 

In section|2j2]we have presented two different approaches to estab- 
lish the stellar parameters characterizing a star. There is still another 
way, based on SpT-Jeff calibrations. 

It is important to note that while the spectral classification of a 
star provides a qualitative description of the stellar characteristics, 
and does not depend of any model calculation, the stellar param- 
eters determination is dependent on the physics used in the stellar 
atmosphere modeling. For example, in the last years there has been 
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an important revision of the cali bration s of stellar parameters of O- 
type stars (see e.g. Martins et al. l2005h . The inclusion of non-LTE, 
line blanketing and wind effects in the stellar atmosphere code cal- 
culations resulted in SpT - T^ff calibrations indicating a lower ef- 
fective temperature for a given spectral type compared to previous 
calibrations (Vacca et al. 119960 . 

The above results (Sects. I2.4I > alerts us that one must use the 
SpT-Teff calibrations with care, because these are model dependent. 
Obviously, this does not mean that calibrations cannot be used at 
all, but it is important to understand that applying a calibration ob- 
tained with one code to create a grid of stellar SEDs for O-type stars 
with a different code would be an inconsistent way of doing. More- 
over, even using the same code, it should be reminded that a certain 
dispersion is always associa ted wi th a given calibration . We r efer 
the reader to Martins et al. boost) and Mokiem et al. j2004h for 
a detailed discussion on the dependence of the SpT - T^ff calibra- 
tions on model assumptions such as microturbulence, metallicity, 
and at mos pheric and wind parameters. As commented by Martins 
et al. ( I200S1 . their theoretical T^g scale should be taken as indica- 
tive since it is expected to have an uncertainty of ±1000 to 2000 
K due to both a natural dispersion and uncertainties inherent to the 
methodology they applied. 



3 PREDICTED EFFECTS OF DIFFERENT SEDS ON 
SURROUNDING NEBULAE 

In this section, we investigate how the differences in the SEDs may 
affect the overall properties of surrounding ionized nebulae. For il- 
lustrative purposes, we have constructed a set of ab initio photoion- 
ization models with cloud-JI, using as an input the SEDs from the 
stellar atmosphere models discussed above. 

To isolate the dependence of the nebular ionization structure 
on the stellar SED from the other factors that may also affect it, such 
as nebular gas distribution and ionization parameter, we limited 
ourselves to spherical constant density models (Nn = 1000 cm"^^, 
Rm =10'^ cm), and the same QQl ") was considered for each set of 
stellar parameters. The adopted nebular abundances were those of 
the Orion nebula and all the models are dust-free. 

3.1 Effect on the nebular temperature 

3.1.1 General considerations 

At each pointin the nebula, the electron temperature, T^, is the re- 
sult of balance between heating and cooling processes. The en- 
ergy gains in a photoionized nebula are usually dominated by 
photoionization of hydrogen atoms, with some contribution of 
photoionization of helium (Osterbrock 1989; see also review by 
Stasinska liooi) . It can be shown that, when ionization equilibrium 
is achieved, the energy gains due to ionization of H " and He " can, 
in a first approximation, be written as: 

G = (< -^^H > ffH+ < £hc > ffHe t), (2) 

"H "He «H 

where < En > is the average energy gained per photoionization of 
an H " atom: 

< En >= J'' — — (3) 

^ We used version 07.02 of Cloudy, last described by Ferland et al. jl998l) . 



Table 5. Nebular temperature averages (over volume) weiglited by the elec- 
tron density predicted by cloudy models using the SEDs from the stellar 
models considered in Sect. 12.31 
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with being the intensity of the ionizing radiation. An expression 
similar to Eq.|3]can be written for < CHe ^■ 

Therefore, < E > depends on the shape of the ionizing SED, 
but not on the stellar luminosity, nor on distance to the ionizing 
source. 

On the other hand, the most important source of cooling in the 
nebula is generally collisional excitation of low-lying energy levels 
of abundant elements (O, S, N) followed by radiative de-excitation. 

3.1.2 Quantitative effects of different SEDs on Te 

Table [5] summarizes various nebular temperature averages (over 
volume and weighted by the electron density) predicted by cloudy 
models using the SEDs from the stellar models considered in Sect. 
lU 

The first thing that can be noted is that, in spite of the large 
differences found between some of the SEDs (see Fig. [T), the re- 
sulting nebular electron temperatures do not differ by more than 
300-350 K. 

To understand this one should note that in the calculation of 
< E > from Eq.[3] the shape of the stellar SED is modified by the 
effect of absorption of stellar radiation in intervening nebular lay- 
ers, e"'"', and by the absorption cross section, ay(ll°). Therefore, 
the contribution of different spectral regions from the SED to the 
nebular heating is not equally important. The functions to be inte- 
grated in the numerator of Eq.|3] for the case of an optical thickness 
of zercQ, are plotted in Fig.|4] with the left panels corresponding to 
H and the right panels corresponding to He. fastwind stellar SEDs 
for stars with Teff = 30000, 35000 and 40000 K (from top to bottom, 
respectively) were used as illustrative. These figures show that the 
average energy gained per photoionization of H" in the nebula is 
independent of the SED above 25 eV for the stars with T^ff = 30000 
and 35000 K, and 35 eV for the T^ff = 40000 K case. Consequently, 
if the SEDs predicted by the stellar atmosphere codes differ be- 
low these energy limits, the nebular electronic temperature may be 

* Inside the nebula, those functions ai'e affected by absorption, and by the 
contribution of the diffuse ionizing radiation produced in the nebula, both 
of which are dependent on the SED only to second order 
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Figure 4. Numerator functions to be integrated to calculate < £h > (left) 
and < Ehc >(right) from Eq.|3] Note the different scales used in each plot. 



different. On the other hand, differences in the SEDs above these 
limits will not have any effect on the heating of the nebular gas. A 
similar argument can be applied to the heating by photoionization 
of He", which is dominated by photons with energies around 30 
eV. Note, however, that the latter is always much smaller than H " 
heating, because of the hhc/wh factor in Eq.|2] 

One could then expect that, for a photoionization model 
with the same nebular abundances, the harder the ionizing fluxes 
for energies between 15-30 eV, the higher the nebular electron 
temperature. However, from inspection of the first column of Table 
[S] it can be seen that actually this is not the general behavior. For 
example, in the Teff = 35000 K case, while cmfgen is resulting in 
the hardest ionizing flux (see Fig.[TJ, the corresponding photoion- 
ization model is the one giving the lowest T^. On the other hand, 
for the same star case, in spite of the fact that cmfgen and tlusty 
are producing a similar heating, there is a difference in between 
the associated cloudy models of 240 K. 

The reason for this is that, although the nebular abundances 
are the same, the various SEDs are producing a different nebular 
ionization structure and hence the ionic fractions of O O S 
S^^, S^, and (main coolants) are different. For example, in 
the first case mentioned above, since the number of O ^ ionizing 
photons is larger in cmfgen (see Fig.[TJ, the population of O will 
be larger in the nebula. Since [O iii] optical and infrared lines are 
very efficient coolants, the overall cooling will be more effective. 
Note that this effect is so strong that the resulting nebular in the 
O region (and hence in the global nebula) is lower than for the 
nebula ionized by fastwind model, which has a weaker ionizing 
flux, but the amount of O ^ ionizing photons is smaller. 

Because of this effect, the shape of the SED in the spectral 
range containing the O ^ and S edges is crucial. This range is 
precisely where the larger dispersion between the various SEDs is 
found, even when the fluxes are compared in the observational ap- 
proach (see Fig.[3](- 

Obviously, the magnitude of these effects depend on the 
star and on the nebula under consideration. For example, in the 
Teff = 30000 K case the O region is very small, and hence the 
global effect of cooling by [Om] lines is not very important. In 
this case, both re(O^) and are similar, and the nebula with a 
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Figure 5. Nebular diagnostic line ratios from the ab initio photoioniza- 
tion models. Models using SEDs from cmfgen, tlusty, v/M-basic, and fast- 
wind are plotted as black, red, blue, and green curves, respectively. Vertical 
lines show the effect a variation of lOOOK in the electron temperature on the 
line ratios. 



higher is the one ionized by tlusty, whose stellar flux is the 
hardest in the 20-25 eV region. 

Figure [5] illustrates the effect of the various SEDs on the 
nebular temperature from a different point of view. Here we plot 
the values of the indicators [S iii] 6312/[S iii] 9532 (left) and 
[N It] 5755/[N ii] 6584 (right) as would be measured through a small 
square slit as a function of the projected distance to the central star, 
i. e. after integrating the line emissivities along a chord through the 
spherical H ii region. Following the same set of colors used along 
the paper, models using SEDs from cmfgen, tlusty, v/u-basic, and 
FASTWIND are represented in black, red, blue, and green lines, re- 
spectively. The electron temperatures that would be derived from a 
given Te-indicator differ by about 200-300 K at most. 



3.2 Effect on the nebular ionization structure 

Contrary to the energy gains, which result from an integration over 
energy of the ionizing photons, the nebular ionization structure de- 
pends on the details of the SED. One can thus suspect that differ- 
ences in the SEDs obtained with various stellar atmosphere codes 
will translate into significant differences in the ionization structures 
of the nebulae. 



3.2.1 General considerations 

The state of ionization of a nebula depends only, in first approxi- 
mation, on the ionization parameter, U, and on the h ardness of the 
ionizing radiation. Following Vilchez & Pagel ( Il988h , 
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Table 6. Integrated ionic ratios (in logarithm) predicted by cloudy pho- 
toionization models using the SEDs from the stellar models considered in 
Sect. [231 



He+/H+ 02+/0+ Ne2+/Ne+ S^*/S^+ Ar^+/Ar* 



are resulting in enormous differences in these values for these quan- 
tities (up to ~ 1 dex for some cases). Obviously, there is a direct 
relation between how these quantities compare and the discrepan- 
cies found in Fig.[T]for the stellar ionizing SEDs. For example, the 
largest difference in the ratio 2(He ") / Q(il°) were found between 
FASTWiND and cmfgen for the 35000 K star case (~0.4 dex) and, 
consequently, the largest difference in the ratio He */H * predicted 
by the photoionization models is also found for the same case (0.35 
dex difference). 

The fact that different SEDs produce different ionization struc- 
tures has important consequences on abundance determinations and 
on the evaluation of the mean effective temperature of the ionizing 
stars. 

Effects on ionic ratios. In Table [7] we consider the nebular ra- 
tios N+/0^ and Ne'^/O"^. A common way to derive N abun- 
dances in Hii regions is to use [Nii] and [On] lines, and correct 
for the unseen N-+ assuming that: N/0 = N+/0^. Similarly, the 
Ne abundance is obtained assuming that Ne/ O = Ne O Both 
formulae are recipes based on similarities of ionization potentials. 
TableQshows that not only these formulae are not necessarily cor- 
rect, but also that the value of these ratios depend on the considered 
SEDs, even the ionization potentials of N + and O ^ are so close to 
each other (as seen in Fig.[TJ. 

Estimation of the mean temperature of the ionizing stars. This 
can be done by considering line ratios arising from different ions 
of the same element. The most popular indicator is the radiation 
softness parameter ;7 = (0^/0^^)/(S^/S^^) proposed by Vilchez 
& Pagel (1988). A better way to estimate the hardness of the 
ionizing radiation field is to plot the observational points in the 
j/(S-Ne)=([S iv]/[S iii ])/([N e iii]/[Ne ii]) vs. [Neiii]/[Neii] plane, as 
shown by Morisset l l2004h . <Tcff> is then derived, together with 
the ionization parameter U, by comparison with grids of models. 
In Fig.|6l we show t he posi tion of the compact Hii region G29.96- 
0.02 (Morisset et al. l2002h with respect to small grids of models 
constructed with each of the four stellar atmosphere codes under 
study. In this case we have excluded the 30000 K star case and con- 
sider additional stellar model^ with Tej =45000 K and log^ = 4.0 
dex. We see that, while cmfgen and tlusty would infer <T^ff> ^ 
40000K, WM-basic would imply <Tsff> ~ 37000K, and fastwind 
would give <T^ff> ~ 38500K 

The differences between the values of <T^ff> obtained using 
different codes are due not only to the difference in the general 
slopes of the SEDs, but also to local properties of the SEDs close 
to the ionization potentials of the involved ions. 

If diagrams like those of Fig. |6] are used to simply order the 
values of <rj.ff> of a sample of Hii regions, one will probably ob- 
tain roughly the same ordering whatever atmosphere code is used 
(provided that one uses a diagram of appropriate metallicity for 
each object<reff>, as described by Morisset 120041) . But, as Fig. |6] 
shows, taking as granted the absolute values of <Teff> is risky, since 
it depends so much on the stellar atmospheres used to determine 
them. 

An additional problem is that the results depend strongly on 
whether the diagrams are constructed using models of dwarf stars 
or of supergiants. In general, for a given T^ff the ionizing SED is 
harder in a star with lower gravity; therefore, lower stellar temper- 
atures will be derived if based on models of supergiants. This is 
the reason why we obtained a larger <reff> for the ionizing star 



' These models were computed following the same ideas as those presented 
in Sect.im 
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Table 7. Results for the N O + and Ne O ^+ ionic ratios predicted 
by CLOUDY photoionization models using the SEDs from the stellar models 
considered in Sect. 12. 31 



(N+/0+)/(N/0) (Ne2+/02+)/(Ne/0) 



CI 1.002 0.086 

Tl 1.000 0.066 

Wl 0.998 0.127 

Fl 1.002 0.139 

C2 0.804 0.085 

T2 0.748 0.132 

W2 0.899 0.157 

F2 0.910 0.192 

C3 0.668 0.385 

T3 0.445 0.345 

W3 0.612 0.767 

F3 0.422 0.317 



"(X-"') Ln^dv Q(X') 

oc u cx U --, (4) 

n(X') C.f^rfv QCH")' ^' 

where n(X ') is the number density of the (-times ionized atoms of 
element X, and Hy is the stellar Eddington flux. Therefore, a ratio 
n(X'^')/;2(X') is, to a first order and for a given U, proportional to 
the relative number of photons able to ionize X', as compared to 
that of Lyman continuum photons. Note, however, that photoion- 
ization models show that Eq.|4]is not exactly fulfilled (Stasiiiska & 
Schaerer 1997) rendering the interpretation more complex. Indeed, 
Eq-Hignores the effect of the modification of the radiation field in 
the nebula due to nebular emission and absorption. 

3.2.2 Quantitative effects on the nebular ionization structure 

Table [6] summarizes the predictions resulting from cloudy models 
using the stellar SED from the various stellar atmosphere codes in 
terms of the ionic fractions (in logarithm) for several elements. A 
comparison of these quantities can give us an idea of the effect that 
differences in the stellar SEDs are producing on the nebular ioniza- 
tion structure. As expected, in some cases photoionization models 



10 S. Simon-Diaz & G. Stasinska 



CMFGEN 



TLUSTY 



1.0 
0.5 
0.0 
-0.5 
-1.0 

-1.5 
-2.0 




1.0 
0.5 
0.0 
-0.5 
-1.0 

-1.5 
-2.0 




-3 -2-10 1 
log [Nelll]/[Nell] (IR) 

WM-basic 



-5 -2-10 1 
log [Nelll]/[Nell] (IR) 

FASTWIND 





1.0 : 


ST 


0.5 r 




0.0 r 








1 


-0.5 r 






-1.0 r 


cn 




O 


-1.5 r 
-2.0 L 





I 

00 



I -4-5-2-1012 
log [Nelll]/[Nell] (IR) 

Figure 6. Diagnostic diagrams proposed by Morisset i2004i) to derive <T^tf> and U by comparing nebular infrared data with results from grids of photoion- 
ization models. The corresponding IR lines used for the diagrams are [S i v] 10. 5/im, [Sin] 18.7/jm, [Neiii] 15.5//m, and [Nen] 12.8//m. The figure shows the 
position of the compact H n region G29.96-0.02 studied by Morisset et al. j2002l) with respect to small grids of models constructed with each of the four stellar 
atmosphere codes under study. The grids were constructed using photoionization models with log U = -3,-2,-1 using fluxes from stellar models with Teff = 35, 
40, and 45 kK (circles, tiiangles and squares, respectively). 
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of G29.96-0.02 than Morisset et al. j2002h . since they used stellar 
models with lower gravities. 



4 SUMMARY 

The aim of this paper was to study how various modem stellar at- 
mosphere codes compare in terms of ionizing SEDs and how this 
affects the emission line spectra of surrounding nebulae. This is im- 
portant, since different codes predict quite different ionizing SEDs, 
which, in turn, impinges on the quantitative interpretation of H ii re- 
gions spectra in terms of formation history, temperature of ionizing 
radiation or chemical composition. 

We have computed models using four state-of-the-art stellar 
atmosphere codes: tlusty, cmfgen, WM-basic, and fastwind. All 
these codes account for non-LTE and line blanketing effects, but 
with different emphasis on specific aspects, tlusty is a plane paral- 
lel, hydrostatic code, while the other three belong to the family of 
models atmosphere codes so-called "unified models". On the other 
hand, while tlusty and cmfgen consider a fully consistent treat- 
ment of line blanketing, wu-basic and fastwind use a fast numeri- 
cal method to account for the line blanketing and blocking effects. 
Three sets of stellar parameters representing a late 0-type (09.5 V), 
a mid 0-type (07V), and an early O-type (05.5 V) dwarfs were 
considered, with the chemical composition and microturbulence 
being exactly the same in all models, so as to perfectly isolate the 
effects on physics and numerical prescriptions from the effect of 
secondary input parameters. 

We have shown that, in order to compare the ionizing SEDs 
predicted by various stellar atmosphere codes, it is important to 



distinguish between what we called evolutionary and observational 
approaches. The evolutionary approach compares SEDs of stars de- 
fined by the same values of T^ff and log g, and is relevant for stud- 
ies based on stellar population synthesis. On the other hand, it is 
the observational approach which is needed when using Hii re- 
gions to probe theoretical SEDs. In this approach, the models to be 
compared do not necessarily have the same T^g and log g, but they 
produce similar intensities and profiles for the stellar H and He i-ii 
lines normally used to estimate the stellar parameters. Note that a 
combined tailored modelling of a massive star and its surrounding 
ionized nebula would precisely use the observational approach. 

We have shown that the study of the optical H and Hei-ii 
stellar lines helps understanding the origin of differences in the 
ionizing SEDs predicted by various stellar atmosphere codes. 
There is a better agreement - in terms of the total number of 
H ionizing photons gCH"), of the ratio Q{Yi&°)IQ{Yl\ and of 
the shape of the SED in the spectral range 13 - 30 eV - between 
models resulting in similar H and He i-ii lines than between models 
constructed with exactly the same T^g and logg. Therefore, to 
probe the ionizing radiation from massive stars by using nebular 
lines from surrounding Hn region!^, it is important to secure high 
quality spectra of the ionizing stars, allowing one to accurately 
determine the stellar parameters. Photoionization modeling will 
then permit one to investigate the ionizing SED above 30 eV 
(where the shape of the SED is known to be very sensitive to the 



^ Table lBll in Appendix|B](accessible on-line only) indicates the available 
optical and IR nebular line ratios useful to test the predicted emergent ion- 
izing SEDs. 
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treatment of line blocking, line blanketing and stellar wind effects). 

We then examined the emission line spectra of ab initio 
models of Hii regions constructed with the various atmosphere 
codes. We found that the nebular temperatures do not differ 
by more than 300-350 K, and that the changes are not always 
in the direction expected from the change in the heating rate: 
differences in the cooling rate, induced by the modification of 
the ionization structure, also play a role. The ionization struc- 
ture itself may change substantially. This has an impact on the 
abundance ratios derived from Hii regions spectra, either by 
tailored photoionization modelling or with the help of ionization 
correction factors. For example, we have found that the values 
of NV 0+ and Ne'V 0^+, commonly identified with the N/0 
and Ne/0 abundance ratios, may differ by 0.2-0.3 dex from 
one SED to another. The changes in the ionization structure 
also affects the temperature of the ionizing stars estimated from 
the emission line ratios produced by the nebula. In the example 
we considered, the estimated T^ff varies from =^ 37000K to 40000K. 

Note that, in this paper, we presented models for dwarf stars 
only. For supergiant stars, the comparative behaviour of models ob- 
tained with the various codes may be different. For example, while 
Morisset et al. (12004) found that the difi'erence between WM-basic 
and CMFGEN models for supergiants with T^g = 40000 K is positive 
in the >35 eV range, we found the opposi te for dwarfs with the 
same T^ff (see also Fig. 8 in Morisset et al. l2004h . It is likely that 
the behavior will change also with the chemical composition. 

The problem is therefore complex and has to be analyzed fur- 
ther. In any case, the limited study presented in this paper shows 
that, than even with the latest generation of stellar atmosphere 
codes, one must be careful when combining them with photoioniza- 
tion codes to establish the ionization correction factors to be used 
in the abundance determination of H ii regions, or to characterize 
the ionizing stellar population from nebular line ratios. 

The follow-up of this study will be to analyze in detail a few 
Galactic H ii regions and their central stars, performing a combined 
stellar and nebular modelling, in order check which of the stellar 
atmosphere codes provides the best ionizing SED. Hopefully, the 
present study and our forthcoming papers will have an impact on 
our understanding of the physics of stellar atmospheres, as well as 
on the reliability of the predictions from stellar synthesis codes and 
of abundance determinations in H ii regions. 
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APPENDIX A: GENERAL CHARACTERISTICS OF THE 
STELLAR ATMOSPHERE CODES 

Below we summarize the main characteristics of the four stellar 
atmosphere codes consider ed in this study (we also refer to Table 
1 in Herrero |2007| and Puls l2008h and give some notes on how the 
specific calculations of the models used here were performed. 

Generally speaking, the four codes account for the nlte and 
line blanketing effects (although they also offer the possibility of 
calculating lte and/or unblanketed models), tlusty is a plane par- 
allel, hydrostatic code, while the other three (fastwind, wu-basic, 
and cmfgen) belong to the family of model atmosphere codes so- 
called "unified models" (Gabler et al. 1989), in which a spherically 
symmetric geometry, with a smooth transition between the photo- 
sphere and the wind is considered. 

Al TLUSTY 

TLUSTY is a plane parallel hydrostatic code widely used nowadays. 
This code is aimed at the study of spectral lines formed in the 
photosphere, where velocities are small and the geometrical exten- 
sion is negligible. The model atmosphere code is completed with 
SYNSPECT, a program for calculating the spectrum emerging from a 
given model atmosphere. 

The basic concepts, equations and numerical methods used are 
described in Hubeny tl988i) . New developments included in subse- 
quent versions are descri bed in detail in Hubeny & Lanz l l 19921 
1 199 A and Hubeny et al. ('1994'). One of the most important fea- 
tures of the new version of the program (v. 200) is that it allows 
for a fully consistent, nlte treatment of metal line blanketing, ei- 
ther by means of Opacity Distribution Functions (ODE) or Opacity 
Sampling (OS). 

The program is fully data oriented as far as the choice of 



atomic species, ions, energy levels, transitions and opacity sources 
is concerned (there are no default opacities built in). Note that to 
calculate realistic line blanketed models with tlusty all the above 
mentioned atomic data need to be included explicitly in the calcu- 
lations. This makes the model calculation to take about 4 hours. 
In order to construct re alistic nlte model atmospheres for O-type 
stars, Lanz & Hubeny l l2003h included about 100000 individual 
atomic levels of over 40 ions of H, He, C, N, O, Ne, Si, P, S, Fe 
and Ni in the calculation of their Ostar2002 grid. The levels are 
grouped into about 900 nlte superlevels. A total of 8000 lines of 
the light elements and about 2 million lines of Fe iii-Fe vi and Ni iii- 
Ni VI are accounted for in the calculations. 

As commented in Sect. 12.11 we used three models from the 
Ostar2002 grid for our study. 

TLUSTY and synspect codes, the Ostar2002 grid and all 
the atomic data files used in it can be downloaded from 
http://nova.astro.umd.edu/ There is also a very complete 
user's guide available. 

A2 CMFGEN 

CMFGEN was originally designed for the analysis of W-R stars and 
LBVs, with very dense winds. Subsequent developments of the 
code have allowed to also extend its application to O stars. Nowa- 
days, CMFGEN can be used for detailed spectroscopic studies of hot 
stars from the (e)uv to the (f)ir ranges. 

The approach driving this code is the inclusion of as many 
lines as possible in order to fully describe the effects of line blan- 
keting, while minimizing the number of level populations to be 
explicitly solve d (using the idea of superlevels, first pioneered by 
Anderson [T989h . Radiative transfer is treated "exactly" in this code 
(all lines are treated in comoving frame, even those from iron group 
elements), meaning that no opacity redistribution or sampling tech- 
niques are used. Consequently, cmfgen is very time consuming 
(model calculation takes anything from 3-12 hours, depending on 
the star and how close the initial guess is from the final solution). 

Since cmfgen was originally designed for the analysis of stars 
with very dense winds, the treatment of the photospheric density 
stratification is considered in an approximated way. At present, the 
velocity and density structures are not self-consistently computed 
but have to be given as inp ut, either as input data or parameterized 
(see in Bouret et al. l2003l an example of how these structures are 
included in the code). 

The CMFGEN models used in this study were kindly generated 
by F. Najarro and F. Martins. They used as initial guess for the 
model calculation previously calculated models with stellar param- 
eters close to the required ones (see Table[T]l. Each model calcula- 
tion took approximately 1 day. 



A3 WM-basic 

WM-basic aims mainly at the prediction of euv/uv fluxes and pro- 
files for hot luminous stars with homogeneous, stationary, and 
spherically symmetric radiation driven winds. In terms of the emer- 
gent spectrum calculation, in -WM-basic models the bound-bound ra- 
diative rates are calculated in Sobolev approximation and the code 
does not include Stark broadening. Therefore, in WM-basic, the op- 
tical lines usually considered for the stellar parameter determina- 
tion are not reliable for diagnostic purposes. 

WM-basic provides a realistic stratification of the density and 
velocity, particularly in the transonic region, through the solution 
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of the hydrodynamics produced by the scattering and absorption of 
Doppler shifted metal hnes. The line contribution to the radiative 
acceleration is parameterized following the formulation by Castor , 
Abbot & Klein ( 1975) — CAK — (see also Pauldrach et al l 19861) 
that uses the concept of line-force multiplier and the force multi- 
plier parameters (K, a and S). 

The authors have developed a fast numerical method which 
accounts for the blocking and blanketing influence of all metal lines 
in the entire sub- and supersonically expand ing atmosphere in a 
consistent manner (see Pauldrach et al. bOOlh . This method makes 
the computational time to be significantly reduced in comparation 
with e.g. CMFGEN calculations. To this end, detailed atomic models 
are implicitly considered by the code for all the important ions. 
These include 149 ionizations stages of 26 elements (H to Zn, apart 
from Li, Be, B and Sc), resulting in a total of 5000 levels. More than 
30000 bound-bound transitions for the nlte calculations and more 
than 4 million lines for the line-force and blocking calculations are 
accounted for. 

Our WM-basic models were generated with the easy-to-use 
Windows interface created by Adi Pauldrach for running the cod^ 
The initial guess of the force multiplier parameters (K, a and 6) 
provided by wu-basic was slightly modified to produce the same 
M and Voo values than the ones used for fastwind and cmfgen mod- 
els (see Table[T]l. The models took ~ 30-50 min. to be converged. 

A4 FASTWIND 

FASTWIND (Santolaya-Rey et al. 1 19971 : Puis et al. 1 20051) is optimized 
for the analysis of optical and IR spectra of OBA-type stars of all 
luminosity classes and wind strengths. 

Among the main objectives that have guided the developers of 
FASTWIND since the initial version of the code a fast performance is 
the highest priority motivation. Since the first version, presented in 
Santolaya-Rey et al. ( Il997h to the latest one (last described by Puis 
et al. 1 2005) . many improvements have been included in the code 
(the main ones are an approximate treatment of the line blanket- 
ing/blocking effects and a consistent calculation of the temperature 
structure by means of the thermal balance of electrons). 

The required computational efficiency is obtained by apply- 
ing appropriate physical approximations to processes where high 
accuracy is not required (regarding the objective of the analysis: 
optical/IR lines). One of these approximations concerns the treat- 
me nt of t he opacity from metals (line blocking). We refer to Puis et 
al. j2005l) for a detailed description of the method. 

The code distinguishes between explicit elements, that are in- 
cluded in detail and can be synthesized later when solving the for- 
mal solution, and background elements producing the line blanket- 
ing and whose occupation numbers are calculated by solving the 
NLTE bound-bound rate equations in Sobolev approximation, in a 
similar way as wu-basic does. The calculation of the line blanket- 
ing effects is hence independent of the elements considered explic- 
itly 

Data from the background ions are taken from Pauldrach et 
al. (l200l|; WM-hasic) and are provided in a fixed form. 

Although FASTWIND Is uot aimed at predicting accurate ionizing 
spectral energy distributions (due to the simplified treatment of line 
blocking), we included it in our comparisons to test the reliability 

' This is distributed as freeware and can be downloaded from the website: 
www.usm.uni-muenchen.de/people/adi/adi .html . 



of the approximate SED predicted by this code, fastwind models 
were calculated with version v8.5 of the code. Since we were only 
interested in the SED and the optical H and He lines, we considered 
H and He as the only explicit elements. This way we could generate 
FASTWIND models in only ~ 10-20 min. 



APPENDIX B: NEBULAR LINE RATIOS BETTER SUITED 
TO TEST THE SHAPE OF THE STELLAR IONIZING SED 

Hii regions emit numerous H and He recombination lines in the 
optical, along with some forbidden lines from abundant metals. 
In addition, a rich spectrum of H recombination lines and fine 
structure lines is observed from the near- to far- IR (Ne, S, Ar: 
-3-40 yum; C, N, O >40 yum). We refer to Stasinska ( i2007l) . who 
gives a list of useful lines from C, N, O, Ne, S, CI, and Ar ions and 
shows how to deal with them. Using appropriate nebular line ratios 
and taking into account the ideas presented above one can impose 
constraints to test the predicted emergent ionizing SEDs. 

Nebular line ratios involving a He i recombination line relative 
to a H I line (e.g. He i A 5876 / H/?, in the optical) are important con- 
straints. In an ionization bounded nebula, such a line ratio depends 
almost exclusively on the ratio 2(He'')/2(H''), provided that he- 
lium is not ionized in the entire nebula (which occurs if the effective 
temperature of the stellar ionizing source is > 40000 K). 

Intensity ratios of lines from successive ions (X'*' and X') of 
the same element can be used to constrain the ratio of ionic abun- 
dances and, as commented in Sect. |3.2.T| the number of stellar pho- 
tons able to ionize X ' relative to the number of Lyman photons, 
GCXO/GCH"). In order to isolate the effect of the SEDs on the ion- 
ization structure from the effect on nebular temperature, it is better 
to use pairs of lines which have a similar dependence on the elec- 
tron temperature, i.e. either pairs of recombination lines, or pairs of 
far IR lines, or pairs of lines with similar excitation potentials (e. g. 
[S III] 9069/[S II] 6716-H30 rather than [S iii] 6312/[S ii] 6716+30). 

Table |BT] gives the list of available nebular line ratios. For N, 
O, Ne, S, CI, and Ar, only nebular and fine structure lines resulting 
from the lower energy levels are considered. In those cases when 
various lines are available in the optical or IR ranges, the most in- 
tense ones are selected, though other possible lines are also indi- 
cated in the last column. 

This paper has been typeset from a TgX/ MTgX file prepared by the 
author. 



14 S. Simon-Diaz & G. Stasinska 

Table Bl. Summary of nebular line ratios useful to test the predicted emergent ionizing SEDs. 



loniz. edge 


Ei„,m (cV) 


^(A) 


Ionic ratio 


Line ratio 


Optical range (A) 


IR range (pm) 


Other lines 


s + 


23.3 


531 


X(S2+)/X(S+) 


[SiiiJ/[SiiJ 


[Sin] 953I/[S u] 6716+30 




[Sin] 9069, 6312 


a* 


23.8 


520 


X(C1 ^+)/X(Cl + ) 


[CI III] /[Chi] 


[Cliii]5517+37/[Clii] 8579 




[Clii] 9123, 6161 


He" 


24.6 


504 


He + /H + 


Hci/Hi 


Hei/15876/H/3^ 






Ai + 


27.6 


448 


X(Ar2+)/X(Ar+) 


[Arm] /[Aril] 




[Ariii]9.0/[Arii]7.0 


[Arm] 21.8 


N + 


29.6 


419 


X(n2+)/X(N + ) 


[N III] /[Nil] 




[N in] 57.3 /[Nil] 121.7 




S2+ 


34.8 


356 


X(S'+)/X(S 2+) 


[Siv]/[Siii] 




[Siv] 10.5/[Sni] 18.7 


[Sin] 33.5 


+ 


35.1 


353 


X(02+)/X(0 + ) 


[0 III] /[On] 


[0 in] 5007/ 10 II] 3726+29 




10 m] 4959, 4363, [0 n] 7320+30 




39.6 


312 


X(Cl3+)/X(Cl2+) 


[CI iv]/[Cliii] 


[Cliv] 8048/[Cliii]5517+37 






Ne+ 


40.5 


304 


X(Ne2+)/X(Ne + ) 


[Nciii]/lNc II] 




[Neiii] 15.6/lNeii] 12.8 


[Neiii] 36.0 


Ar2+ 


40.7 


302 


X(Ar3+)/X(Ar^+) 


[Ariv]/[Ariii] 


[Ariv] 47 11 +40/ [Arm] 7751 




[Ar IV] 7170, [Arm] 7135 


S3+ 


47.3 


262 




[Sv]/[Siv] 








N2+ 


47.5 


261 


X(N3+)/X(N2+) 


[N IV] /[N III] 









